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ABSTRACT 

We consider recent observations of the chromospheric network, and argue 
that the bright network grains observed in the Ca II H & K lines are heated 
by an as yet unidentified quasi-steady process. We propose that the heating is 
caused by dissipation of short-period magnetoacoustic waves in magnetic fiux 
tubes (periods less than 100 s). Magnetohydrodynamic (MHD) models of such 
waves are presented. We consider wave generation in the network due to two 
separate processes: (a) by transverse motions at the base of the fiux tube; and (b) 
by the absorption of acoustic waves generated in the ambient medium. We find 
that the former mechanism leads to an efficient heating of the chromosphere by 
slow magnetoacoustic waves propagating along magnetic field lines. This heating 
is produced by shock waves with a horizontal size of a few hundred kilometers. 
In contrast, acoustic waves excited in the ambient medium are converted into 
transverse fast modes that travel rapidly through the fiux tube and do not form 
shocks, unless the acoustic sources are located within 100 km from the tube axis. 
We conclude that the magnetic network may be heated by magnetoacoustic waves 
that are generated in or near the fiux tubes. 

Subject headings: MHD — Sun: chromosphere — Sun: magnetic fields — - Sun: 
oscillations 
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Introduction 



In the chromosphere on the quiet Sun it is use ful to distinguish betwee n the magnetic 
network on the boundary of supergranulation cells (ISimon fc Leightorull964l ). where strong 
magnetic fields are organized in mainly vertical flux tubes, and internetwork regions in the 
cell interiors, where magnetic fields are weaker and dynamically less important. 

The canonical picture of the magnetic network is that it consists of vertical magnetic 
fields clumped into elements or flux tubes that are located in intergranular lanes, have mag- 
netic field strengths in the kilogauss range, and have diamete r s of the order of 100 km or less 



at their footpoints in the photosphere (e.g.. lGaizauskaslll985l : IZwaanlll987l ). These magnetic 



elements can be identified with bright points seen in images taken in the G-band (4305 A) 
High resolution observations show that these flux elements are in a highly dy namical state 
due to buffeting by convective flows on granular and supergranular scales (e.g., iMuUer et al. 



1994; 


Berger et al. 


1995. 


1998; 



ability of new ground-based telescopes at excellent sites and sophisticated image reconstruc- 
tion techniques it is now possible to examine magnetic elements with an improved resolution 
of ab out 0.17 arcsec and investigate their structure and dynamics in unprecedente d detail 



[e.g. , iBerger et al.l l2004l : iRouppe van der Voort et al.l l2005l : iLangangen et al.l 120071 ) . High- 
quahty observations of photospheric magnetic f ields are now also being obtained with th e 
Solar Optical Telescope (SOT) on Hinode (e.g.. iLites et al.l 120071 . l2008l : iRezaei et al.l 120071 ). 
Magneto-convection models have been developed to understand the three-dimensional (3-D) 
structure and evolution of the magnetic field and its interactions with convective flows (e.g.. 
Voder et allbood : ISchaffenberger et al.l[2006[ ). 



The chrom ospheric network i s most clearly seen in filter images taken in the Ca II H 



fc K lines (e.g., iGaizauskad Il985l : iRuttenI 120071 ) and in the Ca II IR triplet (jCauzzi et al. 



20071 ). In H or K line images the network shows up as a collection of "coarse mottles" 
or "network grains" that stand out against the darker background. The network grains 
are continuously bright with intensities that vary slowly in time, in contrast to the "fine 



mott les" or "cell grains" w. 



[e.g.. 



Rutten fc Uitenbroek 



Fig. 3 in 



lich a re located in the cell interiors and are much more dynamic 



19911 ). In space-time diagrams derived from Ca II H spectra (e.g.. 
Cram &: Damelll983l ). the network grains show up as long bright streaks, indicating 
lifetimes of at least 10 min. The Ca II H line profiles of network grains are more symmetrical 
than those from the cell interior, and bot h the B.2v and B.2R features generally appear bright 
( jLites et al.l Il993l : ISheminova et al.l l2005l ) . Superpose d on this constant brightness are low 



frequency oscillations with periods of 5 to 20 min. (ILites et al.l Il993l : iKontogiannis et al. 



20061 
(e.g. 



Tritschler et al. 



Bloomfield et al. 



20071) : these oscillations have also been observed at UV wavelengths 



20061 ). 
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The physical processes that produce the enhanced emission in Ca II network grains have 
not been identified. Are they heated by wave dissipation, and if so, what is the nature of 
these waves? In this paper we focus on the question why the emission from network grains 
is so constant in time. We discuss the observations, and investigate whether waves with 
periods P < 100 s can provide the heating of the grains. 



Hasan et al.l (120051 . hereafter Paper I) considered (magneto) acoustic waves generated at 
the interface of magnetic fiux tubes and the outside medium. They suggested that such 
waves are generated as a result of interactions of fiux tubes with turbulent downfiows in 
intergranular lanes just below the photosphere. Intergranular lanes have widths of order 100 
km, and downfiow velocities of several km s~^, so velocity fiuctuations (5f ~ 1 km s~^ may 
be expected. Such fiuctuations would produce transverse displacements of the fiux tubes 
within the intergranular lanes. Assuming random displacements with amplitude 6x ~ 10 
km, the typical time scale is r ~ 10 s. Such perturbations would produce transverse waves 
with period P ~ 2r that travel upward along the fiux tube, as well as pressure fiuctuations 
on either side of the tube (Paper I). Therefore, from a theoretical perspective it is reasonable 
to expect the existence of short-period magnetoacoustic waves in fiux tubes (P ~ 20 s). 

In this paper we present models of wave propagation in magnetic fiux tubes that are 
relevant to the heating of the network grains. Specifically, we present two-dimensional (2-D) 
magnetohydrodynamic (MHD) simulations of the interactions of short-period waves with 
magnetic fiux sheets (we shall use the terminology fiux sheet and fiux tube interchangeably). 
Although such simulations do not include the Alfven mode, they provide a possible model 
for the waves involved in heating the Ca II grains. This work is a continuation of the study 
in Paper I. 

The plan of the paper is as follow. In §2 we review recent observations of the chromo- 
spheric network, and we argue that the network grains are heated by an as yet unidentified 
quasi-steady process. We propose that the heating is due to short-period waves (P < 100 
s), although there is little evidence for such waves at present. In §3 we describe a 2-D MHD 
model for the propagation of short-period waves in magnetic fiux tubes. In §4 we present 
simulation results for waves generated inside a fiux tube; interaction of waves on different 
flux tubes are also considered. In §5 we consider the interaction of a flux tube with waves 
generated by acoustic sources located outside the flux tube. The flnal §6 discusses and 
summarizes the main results of this investigation. 
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Observations of Ca II Network Grains 



Figure [T] shows an example of a network element observed with the Dutch Open Tele- 
scope (DOT. iRutten et al.ll2005l ). The different panels show filtergrams in the blue contin- 
uum, G band, Ca II H line center, and Ca II H off band. The DOT uses a filter with a 
1.4 A FWHM passband centered on Ca II H (3968 A), which covers the two chromospheric 
emission peaks on either side of line center, but also includes parts of the inner wings of 
the line, which originate in the upper photosphere. An examination of the G-band image 
confirms the picture of a network patch consisting of multiple flux elements. A comparison 
with the Ca II line center image shows that the excess chromospheric emission is localized 
directly above the photospheric flux tubes, although the bright features seen in Ca II are 
more diffuse than those seen in the G-band. Tim e sequences of such im ages show that the 



chromospheric network is continually bright (e.g.. iTritschler et al.ll2007l ). 



Sheminova et al.l (120051 ) combined high-resolution Ca II H and K spectra of network 



elements with flux tube modeling, and found that the excess brightness in the wings of 
the H and K lines compared with the quiet photosphere is primarily due to low density, 
not to mechanical heating. However, both the B.2v and B.2R emission features appear bright, 
indicating sustai ned non-thermal heating at heights where these features are formed (also see 
Lites et al.lll993l ). The Ca II profiles from network elements are nearly symmetric, indicating 
that the heating does n ot involve strong shocks such as those found in internetwork regions 
dCarlsson fc Steinlll997h . 



RuttenI (120061 . 120071 ) presented reviews and a synthesis of recent high- resolution observa- 



tions of the solar chromosphere. In the Ca II H & K lines, the network shows up as a collection 
of "Ca II bright points" or "grains" . He also identifies an exciting new phenomeno n called 
"stra ws" that extend radially outward from network bright points (see Figure 8 in iRutten 
20071 ). These straws are very thin, occur in "hedge rows", and are very short-lived (10-20 s). 
They appear to be closely related to the so-called type-II spicules recently identified in limb 
observations with Hinode/SOT (IDe Pontieu et al.l l2007bl ). The type-II spicules are much 
more dynamic than ordinary (type-I) spicules, are very thin (width ~ 100 km), have life- 
times of 10-150 s, and seem to be rapidly heated to transition region temperatures, sending 
material though the chromosphere at speeds of 50-150 km s^^. De Pontieu and collaborators 
suggest that type-II spicu l es may be due to small-scale reconnection events in the chromo- 
sphere. iDe Pontieu et al.l (l2007d ) show that the spicules exhibit transverse motions with 
velocities of 10 to 25 km s~^, which provides evidence for the existence of Alfven waves in 
the chromosphere. 

In Ca II observations of network elements on the solar disk, the thin straws are su- 



perposed on the bright grains. According to iRutteru (120061 ): "The network bright points 
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are less sharp than and differ in morphology from G-band bright points ... bright stalks 
emanate a few arcseconds from them at this resolution ... diffuse Ca II H core brightness 
spreads as far or farther" (see his Fig. 2). Rutten conjectures that the bright Ca II network 
grains are nothing but a collection of straws seen nearly end-on. However, this is difficult to 
reconcile with the observed steadiness of the emission from network elements at high spatial 
resolution. Since the straws are very short-lived and dynamic, a superposition of many such 
features would be required to produce apparantly stea dy emission, but s u ch high density of 
straws is not consistent with the limb observatio ns by 



De Pontieu et all (12007 



the plasma in the straws is hot (T ~ 10^ K), as [Rutten J2006h suggests, one would expect 



Also, if 



the straws to produce pure emission profiles, but the observed Ca II H & K line profiles from 
network elements are nearly symmetric and always show a central reversal (see Fig. 5 in 
Sheminova et al.ll2005l ). Therefore, we suggest that most of the emission from Ca II grains 
is due to some quasi-steady heating process at heights of 500 km to 1500 km inside the 
magnetic flux tubes, and is not associated with spicule-like activity at larger heights. 

Our interpretation of the Ca II observations is summarized in Figure [21 where we show a 
vertical cross section of a magnetic network element consisting of several discrete flux tubes. 
We suggest that the Ca II network grains are located inside the magnetic flux tubes, and 
give rise to the bulk of the Ca II emission from the network element. The grains are thought 
to be located at heights between 500 km and 1500 km above the photosphere where the 
flux tubes are no longer "thin" compared to the pressure scale height (about 200 km), but 
are still well separated from each other. The Ca II straws (type-II spicules) have widths 
of order 100 km, and are located at larger heights (several Mm) where the widths of the 
flux tubes are much larger than 100 km. Therefore, we suggest the straws are not directly 
associated with the network grains in the low chromosphere. In Figure [2l we assumed that 
th e straws (type-II spicules) are located at the interfaces between the flux tubes, as suggested 
by [van Ballegooijen fc NisensonI (119981 ). 



Observations have shown "magnetic shadows" surro unding network elements where 



the oscillations in the 2 - 4 minute range are suppressed (jJudge et al.|[200ll : [Krijger et al. 



200l[ : [Vecchio et al.[ [20071 ). These shadows have been attributed to the effects of highly 



structured magnetic field s on wave propagation and mode conversion in the chromosphere 



([Mcintosh &: Judge[[200l[ ). Netw ork elements on the quiet Sun exhibit os cillations with peri- 



ods in the range 5 - 20 min. (e.g.. [Lites et al.[[l993[ : [Bloomfield et al. hood ). These long-period 



oscillations are thought to be driven by the interaction of vertical or inclined magnetic fiux 
tubes with global p-modes and convective flows in the photosphere, producing magnetoacous- 



(De Pontieu et al. 


2004 


2005 


Hansteen et al.|[2006: 


Heeeland. De Pontieu & Hansteen 


2007; 


Rouppe van der Voort et al. 


2007 


)• 


Jefferies et al.l ( 


2006 


) argued that inclined magnetic 
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fields at the boundaries of supergranules may act as portals through w hich low-frequency 
magnetoacoustic waves can propagate into the solar chromosphere (also see lMcIntosh fc Jefferies 
2006h . 



It seems unlikely that long-period waves are also responsible for the heating of the Ca II 
network grains in the low chromosphere. Simulations of shock waves with periods P ~ 200 
s in a plane-parallel, non-magnetic atmosphere have shown that such waves produce large 
asymmetries in the Ca II H line profiles, and strong variations in the integrated emission. 
If the grains were heated by such long-period waves, they should exhibit similar strong 
intensity variations. This is not observed, so the long-period waves observed in network 
elements cannot be the main source of heating for the Ca II grains. 

Network grains could possibly be heated by dissipation of waves with shorter periods 
(P < 100 s). Ground-based observations of high-frequency waves in small network elements 
are affected by seeing, so it is possible that waves with periods P < 100 s do exist in 
network elements but are simply not observable fr om the ground. This hypothesis could b e 
tested using Ca II H images from Hinode/SOT (ITsuneta et al.l 120071 : iKosugi et al.l 120071 ). 
keeping in mind of course that the passband of the Ca II H filter also includes a significant 
photospheric contribution. Alternatively, network grains cou ld be heated by dissipat ion of 
Alfven waves. There is evidence for Alfven waves in spicules ( iDe Pontieu et al.ll2007d ). and 
if these waves originate below the photosphere they must travel through the heights where 
the Ca II network grains are located. At present there is no direct observational evidence 
for Alfven waves in flux elements in the photosphere, nor at heights where the Ca II H line 
is formed. Also, simulations of such waves in flux tubes would require 3-D MHD models. 
Therefore, in the present paper we focus on magnetoacoustic waves, which can be simulated 
with 2-D MHD models. 



3. Models for Magnetoacoustic Waves in Network Elements 

3.1. Background 

Earlier work on w ave dynamics idealized the network in terms of thin flux tubes (e.g., 
Roberts fc Webblll978l ). and treated wave prop agation in te rms of the well known transverse 
(kink) and longitudinal (sausage) modes (e.g.. ISpruitlll98ll ). Several investigations have fo- 
cused on the generation and propagat ion of transverse and longitudinal wave rnodes and their 
dissip ation in the chromospher e fe.g..lZhugzhda et al.lll995t iFawzy et al.ll2002l : lUlmschneider 
and references therein). iHasan fc KalkofenI (119991 ) examined the excitation of trans- 



2003 



verse and longitudinal waves in magnetic flux tubes by the impact of fast granules on flux 
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tubes, as observ ed by iMuUer &: Roudierl (119921 ) and iMuller et al.l (jl99J) , and following the 
investigation by Choudhuri et al.r ( 1993 ). who studied the generation of kink waves by foot- 
point motion of flux tubes. The observational signature of the modelled process was highly 
intermittent in radiation emerging in the H and K lines, contrary to observations. By adding 
waves that were generated by high-frequency motion due to the turbulence of the medium 
surrounding flux tubes the energy injection into the gas inside a flux tube became less inter- 
mittent, and the time variation of the emergent radiation w as in better agreem ent with the 
more steady observed intensity from the magnetic network (IHasan et al.ll2000l ). 



The above studies on the magnetic network make use of two important idealizations: 
they assume that the magnetic flux tubes are thin, an approximation that becomes ques- 
tionable at about the height of formation of the emission peaks in the cores of the H and K 
lines; and they ne glect the interact i on of neig hboring flux tubes. S ome progress on the first 
issue was made by [Rosenthal et al.l (120021 ) and lBogdan et al.l (120031 ) . who studied wave prop- 
agation in a 2-D stratified atmosphere, assuming the magnetic field can be approximated by 
a potential field. They examined the propagation of waves that are excited from spatially 
localized sources in the photosphere, and found that there is strong mode coupling between 
fast and slow waves at the so-called magnetic canopy, which they identify as the surface 
where the magnetic and gas pressures are comparable. However, potential-field models do 
not provide an accurate description of the boundary layers between the flux tubes and the 
ambient medium, where electric currents flow perpendicular to the magnetic field. In Pa- 
per I we considered more realistic, magnetostatic equilibrium models in which flux sheets are 
embedded in a field-free medium. In such models the boundary layer between the flux tube 
and its surroundings is sharper than in a potential field model (the width of the boundary 
layer is spatially resolved, see ^3.2). This st rongly affects the propagation and reflection 
of waves at the boundary (IHasan et al.ll2006[ ). Also, the shape of the {3 = 1 surface in a 
magnetostatic model can be quite different from that in a potential-field model. 

In the present paper we further develop the model for MHD wave propagation presented 
in Paper I. The initial configuration consists of one or more fiux sheets in magnetostatic 
equilibrium with a relatively sharp interface between the fiux sheets a nd the surrounding gas 



Rosenthal et al 



2002 



unlike previous stu dies in which a potential field is considered (e.g 
Bogdan et al.ll2003l ). The present work goes beyond Paper I in three important aspects: 
firstly, it seeks to understand how the dynamics of a fiux tube is affected when the magnetic 
field in the wave excitation region is not that strong, so that the j3 = 1 surface (approximately 
the level where the sound and Alfven speeds become equal) lies above the base on the 
tube axis. Observations have revealed that the magnetic field strengths in the network are 
not uniformly kilo gauss but have a range of values that go from strong to moderate (e.g., 
Berger et al.ll2004l ). For moderate field strengths, wave generation occurs in a region where 
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the field is essentially weak. Such waves undergo transformation higher up in the atmosphere. 
Secondly, we allow for the presence of multiple flux tubes in a network patch. We investigate 
the dynamical coupling between neighboring flux tubes and identify the regions responsible 
for the enhanced emission in the magnetic chromosphere. Finally, we examine the hypothesis 
that the magnetic network can be heated by acoustic waves from the ambient medium. More 
details of the model are provided in the next subsections. 



Hansteen et al.l (120061 ) and lDe Pontieu et al.l (j2007al ) performed advanced 2-D numerical 



simulations that span the entire solar atmosphere from the convection zone to the lower 
corona. These simulations describe the evolution of a radiative MHD plasma: the model 
includes the effects of non-gray, non-LTE radiative energy transport in the photosphere and 
chromosphere, as well as field-aligned thermal conduction and optically thin radiative losses 
in the transition region and corona. They find that MHD waves generated by convective 
flows and oscillations in the lower photosphere naturally leak upward into the magnetized 
chromosphere, where they form shocks that can drive spicule-like jets. These models include 
more physics than the one presented here, but we beheve our approach is adequate for the 
task of simulating short-period magnetoacoustic waves. 



3.2. Equilibrium Model 

In the solar photosphere the magnetic flux tubes are spatially separated from each 
other, and are rooted in intergranular lanes. In Paper I we considered a model for the 
structure of a magnetic network element as a collectio n of unipolar flux tubes (also see 



Cranmer fc van Ballegooijenl l2005l : ICranmer et al.l 120071 ) . We presented a two-dimensional 



model with an array of identical flux tubes separated by 1200 km. The flux tubes expand 
with height and merge into a monolithic structure at a height of about 500 km. In the 
present paper we present modified versions of this model. The method for constructing such 
models is summarized below. 

The equilibrium magnetic field B(a;, z) is a function of horizontal coordinate x and height 
z above the level where the optical depth rsooo = 1 in the photosphere. The field is assumed 
to be periodic in the x coordinate; the period length equals the distance between the flux 
tubes. The tubes are assumed to symmetric with respect to their vertical axes, therefore in 
the equilibrium model only one half of one flux tube needs to be considered. For the model 
presented in §4.1, the size of the computational domain is given by Ax = Lx/2 = 1800 km 
and Az = 2950 km, and we use a grid spacing of 5 km. At the lower boundary {z = 0), the 
flux tube has a nearly Gaussian profile with 1/e-width of 65 km. The 2-D structure of the 
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flux tube is determined by solving the force balance equation: 

- Vj9 + pg + -^(V X B) X B = 0, (1) 

where g = —gz. is the gravitational acceleration, z) is the gas pressure, and z) is the 
density. The third term describes the Lorentz force due to electric currents at the boundary 
between the flux tube and its surroundings. The magnetic field is written in terms of a flux 
function A{x, z) such that = ~dA/dz and = dA/dx. The gas pressure is given by 

p{x,z)=pUz)[l + (3,'FiA)], (2) 

where Pint{z) is the internal gas pressure as function of height along the axis of the flux tube 
(y4 = 0); (3o is the ratio of gas- and magnetic pressures at the base (on axis); and F{A) is a 
function describing the variation of gas pressure across field lines (F = on axis and F = 1 
in the ambient medium). A similar expression holds for the gas density. The functions Pint{z) 
and F{A) are described in Paper I, except that in the present paper we use /3o = 2. Equation 
([1]) is solved for A{x,z) using an energy minimization technique (see Paper I). Table 1 lists 
the basic parameters of the equilibrium model on the tube axis and in the ambient medium 
at the base {z = 0) and at height z = 1500 km in the chromosphere. 

To construct the initial conditions for the MHD dynamical model, the equilibrium model 
is mirrored horizontally to create a single flux tube, and then further replicated to produce 
an array of flux tubes. Only the lower part of the equilibrium model is used {z < 1600 km). 
Therefore, the magnetic field at the top of the MHD model is not vertical. 



3.3. Dynamical Model and Boundary Conditions 



We consider wave generation in the configuration described in the previous section by 
perturbing the lower boundary and solving the 2-D MHD equations in conservation form 
for an inviscid adiabatic fluid. These consist of the usu al continuity, morn entum, entropy 
(without sources) and magnetic induction equations (see lSteiner et al.lll994j ). The unknown 
variables are the density, momentum, entropy per unit mass and the magnetic field. We 
assume that the plasma consists of fully ionized hydrogen with a mean molecular weight of 
1.297. The temperature is computed from the specific entropy and the pressure is found 
using the ideal gas law. 



The MHD equations are solved follow i ng th e numerical procedure given by lSteiner et al. 



(I1994J ) and also described in iHasan et al.l (120051 ). Briefly, the equations are discretized on a 
2-D mesh and solved using a finite difference method based on the flux-corrected transport 
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(FCT) scheme of lOran fc Boria (Il987l ). The time integration is exphcit and has second-order 
accuracy in the time step. Small time steps are required to satisfy the Courant condition in 
the upper part of the domain, where the Alfven speed is large. 

Periodic boundary conditions are used at the horizontal boundaries. At the top bound- 
ary, (a) the vertical and horizontal components of momentum are set to zero; (b) the density 
is determined using linear logarithmic extrapolation; (c) the horizontal component of the 
magnetic field and temperature are set equal to the corresponding values at the preceding 
interior point. The vertical component of the magnetic field is determined using the condi- 
tion V • B = 0. Similar conditions are used at the lower boundary, except for the density, 
temperature and horizontal component of the velocity (or momentum). The density is ob- 
tained using cubic spline extrapolation, the temperature is kept constant at its initial value 
and the horizontal or vertical velocity is specified at the lower boundary. 



3.4. Limitations of the Model 

Before describing the simulation results, it is appropriate to mention some of the limita- 
tions of the present investigation. Our model does not include the convection zone, therefore 
the wave excitation mechanism is somewhat idealized. Also, we have used an adiabatic en- 
ergy equation, thereby neglecting radiative losses. The top boundary conditions cause wave 
reflection, so we can run the simulations only for a short period of time (about 180 s). Also, 
with the present code we are limited in our choice for the height of the upper boundary; for 
larger heights the Alfven speed near the top becomes too large and the time step too small 
to make useful simulations. 

Our analysis is based on a 2-D treatment in which the Alfven wave is absent, and the 
coupling of the Alfven mode to the slow and fast modes is neglected. Furthermore, in three 
dimensions the acoustic waves generated by transverse motions of a flux tube can travel 
azimuthally around the flux tube, which is not possible in a 2-D flux sheet. Therefore, the 
pressure fluctuations 5p surrounding a 3-D tube should be smaller than those near a 2-D 
flux sheet. However, based on a simple model for waves generated by an oscillating rigid 
cylinder in a uniform medium, we estimate that the pressure fluctuations are not significantly 
reduced compared to the 3-D case when the wave period P < 2'kR/Cs-, where R is the flux 
tube radius and Cg is the sound speed in the ambient medium. For R = 100 km and = 7.1 
km s~^ we require P < 88 s, so our approximation of the flux tubes as 2-D sheets is adequate 
for short-period waves, but would not be accurate for longer wave periods. In this paper we 
consider waves with a period of 24 s. 
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4. Results for Wave Generation Inside the Flux Sheets 

4.1. Dynamics of a single flux sheet 

We first consider wave generation in a single flux sheet driven by transverse periodic 
motions at the lower boundary. Compared to Paper I we consider a lower value of the 
magnetic fleld strength of 1000 G on the axis at the base z = 0. In this case, the (3 = 1 
surface occurs above the base within the sheet as opposed to Paper 1 where the fleld at the 
lower boundary is 1500 G (on the ax is). The reason for treating this case is the evidence 



from observations (IBerger et al.ll2004l ) for the existence of network elements with moderate 
fleld strengths. We examine the wave pattern that develops in a flux tube due to a transverse 
motion at the lower boundary given by the following form: 

v,ix,0,t) =vosm{2TTt/P), (3) 

where vo denotes the amplitude of the horizontal motion and P the wave period. Similar to 
Paper 1, we choose vq = 750 m s~^ and P = 24 s. Figures [3^ and [St shows the temperature 
perturbation AT (the temperature difference with respect to the initial value at the same 
location) at 75 s, 122 s and 153 s. The black curves denote the magnetic fleld lines and the 
white curves depict constant (3 contours with values 0.1 (upper curve), 1.0 (thick curve) and 
10.0 (lower curve). 

The horizontal motions at the lower boundary produce compressions and decompressions 
of the gas in the flux tube which generate an acoustic like wave (most effectively at the 
interface between the tube and ambient medium as shown in Paper 1) that propagates 
isotropically with an almost constant sound speed. This can be discerned by the almost 
constant spacing in the semicircular color pattern. In the central section of the tube, a 
transverse slow MHD wave is generated that is essentially guided along the fleld lines. At 
t = 75 s, we flnd from Figure [3^1 that the wave pattern is conflned below the /3 = 1 surface 
{z = 0.5 Mm). In this region, where the magnetic fleld can be regarded as weak, the acoustic 
(fast) mode travels ahead of the (slow) MHD wave (which travels at the Alfven speed). On 
the tube axis the acoustic and Alfven speeds at the base are 7.1 km s~^ and 5.6 km s~^ 
whereas at z = 500 km they are 7.3 km s~^ and 8.0 km s~^. At the /3 = 1 level there is 



a strong coupling between the two modes ( as previously demonstrated by [Rosenthal et al. 



2OO2I : iBoedan et al.ll2003l : iHasan et all 120051 ). 



Figures (Hb and ^ depict the situation at instants where the wave has crossed the P = I 
surface. Above this level, the acoustic wave (which is now a slow mode propagating along the 
fleld) continues to travel at about 8 km s~^ but the MHD wave (fast mode) is speeded up due 
to the rapid increase of the Alfven speed with height. The faint blue and yellow halos in the 
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flux tube above the last semi-circle shows the fast mode. It should also be noted that as the 
acoustic-like waves move upwards, their amplitude increases due to the decrease in density 
with height resulting in compressional heating. From Figures [Hb and [3b we find temperature 
enhancements of around 900 K in the upper atmosphere due to amplitude increase of the 
acoustic wave amplitude. 

A comparison of the present calculations for a network flux element with a field on the 
axis of about 1000 G (corresponding to (3 = 2) with the case treated in Paper I of a stronger 
field of 1500 G {(3 = 0.5) shows that there are several qualitative similarities in both cases. 
These include the acoustic pattern in the ambient medium and weak field regions [jS ^ 1) 
as well as the strong heating in the central regions of the tubes. The main difference is in 
the character of the waves near the tube axis: in Paper I, the field in the central regions 
is strong at all heights, whereas in the present case, the field becomes strong only above 
z = 0.5 Mm leading to a layer where strong mode coupling occurs. It is wort h while to 



clarify what happens to a wave when it crosses this layer. As pointed out by ICallyi (120071 ). a 
fast wave in a high (3 region (essentially acoustic in character) propagating along a field line 
changes labels from fast to slow as it crosses the (3 = 1 (or more accurately the one on which 
the sound and Alfven speeds are equal, but this difference is not significant for the purpose of 
the present analysis). In the /3 < 1 region, this wave is now a slow mode, but still acoustic in 
nature. This phenomenon is essentially mode transmission. However, the incident fast wave 
can also be transformed to a magnetic mode (a fast wave in the magnetized region). This 
corresponds to a mode conversion. In practice, when the wave propagation vector makes a 
finite angle with respect to the magnetic field lines, both mode transmission and conversion 
occur. 

Thus, in Paper I the wave generated in the tube due to footpoint motion is a fast 
transverse MHD mode that continues to be fast (near the axis) as it propagates upwards. 
On the other hand, in the present case, the footpoint motions produce a slow transverse 
MHD wave that is transmitted as a fast MHD (still transverse) as it crosses the (3 = 1 layer. 
There is, however, no mode conversion. 



4.2. Dynamics of multiple fiux sheets 

Let us now extend the analysis to a situation where we have multiple fiux sheets. This 
allows us to incorporate the infiuence of neighboring fiux tubes on each other. Figure H] 
shows the wave pattern arising due to a transverse motion of the lower boundary where 
the horizontal component of the velocity has the form given by equation ([3]) with the same 
parameters as before. The fiux tubes in Figure H] are identical to each other and have a 
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field distribution on the base that is the same as in Figure [31 However, above the base the 
field lines do not flare out as much as in the previous case, but straighten out at a lower 
height due to the effect of adjacent flux tubes. At t = 75 s, the waves in the individual 
tubes are sufficiently well separated from each other and the wave pattern in each tube 
is qualitatively similar to the case treated in §4.1. However, at t = 122 s waves emanating 
from neighboring tubes interact with each other especially in the ambient medium. However, 
the wave pattern in any tube is not significantly affected by the presence of its neighbors. 
Furthermore, the slow magnetoacoustic waves above the /? = 1 surface are confined close to 
the central regions of the tubes, where they steepen and produce enhanced heating. This 
heating appears to be dominantly caused by the wave motions generated at the footpoints 
and not by the penetration of acoustic waves from the ambient medium or by waves coming 
from neighboring tubes (we shall examine this aspect in greater detail in §5). A companion 
video of the simulation corresponding to Figures|l](a)-(c) is shown in Videol.mpg. In addition 
to the temperature perturbation AT, we also display velocity vectors (shown by arrows). 
The animation enables one to discern the separation of the modes in the strongly magnetized 
region {(3 <C 1). The transverse MHD wave can be easily identified in the central sections 
of the tubes. The first front of this wave crosses the {3=1 contour at about 70 s and is 
speeded up as it moves into the higher atmosphere, reaching the top boundary at about 105 
s. On the other hand, the slow magnet acoustic wave (for which the velocity is principally 
aligned with the magnetic field) reaches the top boundary at about 180 s. The final frame of 
the video at t = 174 s clearly shows the presence of localized shock-like features with strong 
temperature enhancements of about 3000 K and fiows close to the sound speed. However, 
the temperature perturbations and fiows in the ambient (weak field) medium are very small. 

It is instructive to compare the velocity pattern for the configurations treated in the 
previous two cases in order to identify the nature of the waves in different regions of the 
atmosphere. We decompose the flow into components along the fleld {Vg) and normal to it 
(Ki). Positive (negative) values of Vn correspond to the velocity component parallel (anti- 
parallel) to X-axis when the fleld is vertical. Figures [5^ and [5)3 show the fleld aligned 
{Vs) and transverse (V^) components of the velocity for a single tube at t = 120 s for the 
same parameters as before. We only consider regions with a fleld strength greater than 
20 G (for weak flelds the decomposition into longitudinal and transverse components is not 
meaningful) . 

Figures [St and [5)1 show similar diagrams for the model with multiple tubes. In Figs. [5^ 
and [Sb we see an antisymmetric (about the axis) velocity pattern in the strong fleld region 
{(3 <^ 1); this pattern can be identifled with the longitudinal acoustic mode. On the other 
hand, in Figs. [Sjo and[5li we see a symmetric pattern as a uniform halo, corresponding to the 
transverse MHD mode. There is a clear separation of the modes above the (3 = 1 level, where 
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the fast mode speed increases rapidly with height. The associated increase in wavelength 
manifests itself as an increase in the spacing of the color contours. The shape of the contours 
is a consequence of the Alfven speed being maximum on the axis (at a given height). Figures 
^ and [Hb give a close up of the temperature perturbation AT and velocity pattern near 
the central region of the tube, corresponding to the multiple and single tube configurations, 
respectively. Both cases clearly show the change in the flow pattern close to the /3 = 1 level 
and the anti-symmetric structure of the acoustic waves away from the tube axis. 

We find that the wave patterns in the single and multiple tube cases are qualitatively 
similar. In both cases the acoustic and magnetic type of modes are present. Even though the 
initial perturbation at the base {z = 0) is transverse, strong field aligned motions develop 
particularly close to the central region of the tube, similar to the finding in Paper I for a 
single tube, which suggests that this is a general feature of the network. 



5. Results for Wave Generation in the External Medium 

So far we have examined the excitation of waves in flux tubes due to transverse motions 
of their footpoints. We now investigate the situation when the source of waves is in the 
field free medium, and we consider how the waves interact with the flux tube in order to 
assess how effectively they penetrate into flux tubes and heat the network. Let us consider a 
localised source that generates acoustic waves due to vertical periodic motions. The acoustic 
source is located at the lower boundary of the model {z = 0) and has a Gaussian profile in 
X with a 1/e- width of 50 km. As before, we take the amplitude of the motions as 750 m s~^ 
and period as 24 s. As we shall see the nature of wave excitation in the tube depends on the 
"angle of attack" or the angle between the direction of the wave front and a field line. We 
treat two cases corresponding to distant and near acoustic sources on the lower boundary. 

First consider an acoustic source at a; = 2.3 Mm, which is 500 km away from the axis 
of the flux rope. Figures [7](a)-(c) show the temperature perturbations at different epochs, 
corresponding to t = 76 s, 101 s and 140 s respectively. It is seen that the acoustic wave 
strikes the tube (/? ~ 1 surface) almost normal to the field and generates a fast MHD wave 
in the tube; for a more thoroug h description of what happens at this interface, see the lower 



panel in Fig. 3 of ICallyi (120071 ). There is also slow magnetoacoustic wave generated that 
essentially travels along the field lines. Figures [7](b)-(c) shows the distortion of the wave 
front as the wave propagates into the low (3 region and where it is refracted upwards (away 
from the tube axis). A more informative way to examine the nature of the waves is to plot 
the longitudinal and transverse velocity components, which is shown for the present case 
in Figures [H](a)-(b) at time t = 140 s (for 5 > 20 G). The field aligned component Vg, 
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that essentially depicts the acoustic wave, does not penetrate significantly into the central 
regions of the tube as we see from Figure [8]^a). On the other hand, the normal component 
Vn, associated with the fast MHD mode (where the (3 < 1) dominates in the central regions 
of the tube. This wave is generated due to mode conversion when the acoustic wave strikes 
the P = I surface. 

Now consider an acoustic source at x = 1.9 Mm, only 100 km from the fiux tube axis. 
The vertical field Bz{x,0) at the base has a nearly Gaussian profile with 1/e- width of 65 
km, and the acoustic source has 1/e- width of 50 km, so there is significant overlap between 
the two profiles. Figures [7|^d)-(f) show the temperature perturbation for t = 76 s, 101 s and 
140 s, respectively. In contrast to the previous case, such a source excites waves that travel 
with almost no distortion by the fiux tube, as can be seen by the fact that the wave front 
remains close to circular. From Figures [7](e) we find that the angle between the wave vector 
and the field lines is very small at the /? = 1 level. In this case, there is only a small amount 
of mode conversion into a transverse magnetic mode. In the magnetized region {B > 20 G), 
the waves are dominantly field-aligned, as can be seen in Figures [8](c)-(d) at time t = 140 s, 
which show the field-aligned and normal components of the velocity. In the region with 
P < 1 these waves are slow magnetoacoustic modes, and the fast mode has much less power. 



6. Summary and Discussion 

Recent observations of the chromospheric network were discussed. We argued that Ca II 
network grains represent plasma located at heights between about 0.5 Mm and 1.5 Mm, and 
are distinct from the spicule-like features observed at larger heights (see Fig. [2]). The heating 
in network grains must occur in a sustained (as opposed to intermittent) manner. The physi- 
cal processes involved in this heating ar e not well understood. L ong-period waves are thought 



to be involved in driving spicules (e.g.. lDe Pontieu et al.ll2004l ). but we argued against long- 



period waves as the cause of the heating in Ca II network grains. Instead, our hypothesis 
in this paper is that the heating in network grains is caused by weaker shocks that occur 
repetitively at short time intervals {P < 100 s) to maintain the enhanced temperature of the 
network grains. The weakness of the shocks is needed to explain the observed symmetrical 
shapes of the Ca II H & K line profiles from network grains. The symmetry of the line 
profiles and relative constancy of the Ca II emission may be explained if the shocks have 
small horizontal scales (a few hundred km) and multiple shocks are superposed along the 
line of sight. At present there is no observational evidence for such short-period waves in 
network elements. 



We studied the propagation of short-period magnetoacoustic waves in magnetic ele- 
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ments, using the 2-D MHD model presented in Paper I. The waves are launched as transverse 
kink waves at the base of the photospheric flux tubes. The present work complements Paper I 
in two important respects: (a) we consider an equilibrium model with lower magnetic fleld 
strength, so that the (3 = 1 surface is located in the upper photosphere; (b) we examine the 
interactions of waves from neighboring tubes. We have identifled an efficient mechanism for 
the generation of longitudinal motions in the chromosphere from purely transverse motions 
at the base of tubes. These longitudinal motions are associated with slow magnetoacoustic 
waves that can contribute to a localized heating of the chromosphere. 

The qualitative features of the results for a single flux tube are similar to those found 
in Paper I. Irrespective of whether the fleld at the base of the flux tube is strong or weak, 
the transverse waves at the base are converted at chromospheric heights into two sets of 
slow-mode shocks that are 180° out of phase with each other and propagate upward along 
the magnetic field lines on opposite sides of the fiux tube axis. We propose that network 
grains are heated by dissipation of such shocks. This heating occurs mainly in the central 
part of the flux tube. We also considered an equilibrium conflguration with multiple tubes, 
and found that through the interference of waves emanating from neighboring tubes, there 
is some heating at the tube interfaces. However, the dominant heating occurs near the tube 
axes. 

For the case of a single flux tube we found that the width of the shocked region at time 
t = 153 s is about 1200 km (see Fig. 3c), a signiflcant fraction of the assumed width the flux 
tube at large height. For the case with multiple flux tubes, the width of the shocked region 
is about 900 km, or 3/4 of the width of one flux tube (see Fig. 4c). Therefore, the spatial 
extent of the shocked region in our models is larger than the subarcsecond size of the Ca II H 
network grains observed by DOT (see Fig. Ic). Part of the difference may be ex plained by the 



fact t hat the Ca II H images contain a signiflcant photospheric contribution (jRutteru 12006 



20071 ). which originates at lower heights where the flux tubes are thinner. Another possible 
explanation lies in the details of spectral line formation: as the shocks propagate to larger 
height, the density behind the shocks decreases and the Ca II H line source function drops 
below the Planck function, producing the "absoption" feature at line center. Therefore, the 
chromospheric contribution to the observed emission comes from a limited range of heights, 
and may reflect only the lower part of the shocked region shown in Figures 3c and 4c. To 
determine whether the proposed models can explain the observed widths of network grains, 
more detailed modeling of the formation of the Ca II H line is needed. 

We also investigated the interaction of a flux tube with acoustic waves from the external 
medium. A careful investigation shows that their penetration int o the tube depends on the 



angle between the wave front and fleld lines, or "angle of attack" (ICarlsson fc Bogdanl 12006 
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Callyl 120071 . and references therein) . ISchunker fc Callyl (120061 ) and ICallyl (120071 ) examined 



the interaction between an acoustic wave and a magnetic field, and found that at the level 
where the sound and Alfven speeds are equal there is both mode transmission and mode 
conversion. The former refers only to a nomenclature change as a fast wave in the ambient 
medium is transmitted as a slow wave in the fiux tube. However, the transmitted wave 
continues to be acoustic in character. On the other hand, in the latter case, the wave entering 
the magnetic medium remains fast throughout and changes character from acoustic in the 
ambient medium to magnetic in the fiux tube, which corresponds to a mode conversion. We 
examined this phenomenon by considering two limiting cases of distant and near acoustic 
sources. Our calculations clearly show that acoustic waves, generated in the ambient medium, 
do not significantly penetrate inside the fiux tube unless the acoustic source is located in 
its immediate vicinity (i.e., within about one fiux tube radius). Consequently, short-period 
waves from the external medium do not offer a promising mechanism for heating the magnetic 
network at heights below 1500 km. 

Future observations with Hinode/SOT should be able to establish whether short-period 
magnetoacoustic waves exist in flux elements at the heights where the Ca II B.2v and B.2R 
emission features are formed, and if so, whether they play a signiflcant role in chromospheric 
heating. If such waves are not found, Alfven waves may provide an alternative, but such 
waves may be difficult to detect in solar disk observations. Future modeling should consider 
the interactions of waves with magnetic flux tubes in three dimensions in order to allow for the 
coupling between all three MHD wave modes. The models should include a more realistic 
energy equation than the one considered here. Both p-mode waves and waves generated 
locally at the base of the flux tubes should be considered. Based on such models, Ca II 
H and K proflles should be simulated to determine whether the models can reproduce the 
observed symmetrical proflles from network grains. 



We thank the referee for detailed comments on an earlier draft of this paper. The Dutch 
Open Telescope is operated by Utrecht University at the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astroflsica de Canarias. 
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Table 1. Parameters on the Tube Axis & Ambient Medium 







Tube Axis 




Ambient Medium 


Variable 




z = 1500 km 


z^O 




z = 1500 km 


Temperature 


4700 K 


8700 K 


4700 K 




8700 K 


Density 


2.6 10"^ g cm-3 2.1 10-^^ g cm'^ 


3.9 10- 


7 — S 

g cm 


3.2 10^12 g 


Pressure 


8.0 10^ dyn 


cm~^ 2.2 dyn cm~^ 


1.2 10^ 


dyn cm-2 


3.3 dyn cm-^ 


Sound speed 


7.1 km 


13 km s~^ 


7.1 km 


s-i 


13 km s"^ 


Alfven speed 


5.6 km s~^ 


85 km 








Magnetic field 


1000 G 


44 G 








/3 


2.0 


2.8 10-2 









-23- 





Fig. 1. — Dutch Open Telescope (DOT) observations of a magnetic network patch in the 
solar photosphere and chromosphere. The images were taken at solar disk center on 2007 
April 12 at 9:35 UT. (a) blue continuum; (b) G band at 4305 A; (c) Ca II H hne center; 
(d) Ca II H offband (-2.35 A). The network patch consists of many small flux elements of 
sub-arcsecond size. The bright network grains seen in the Ca II H line center image coincide 
with G-band bright points, indicating that the chromosphcric heating is localized directly 
above the photospheric flux elements. The Ca 11 network grains are more diffuse than the 
G-band bright points, which may be due to the spreading of flux tubes with height. 
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Fig. 2. — Scliematic diagram showing the structure of a magnetic network element on the 
quiet Sun. The thin half-circles at the bottom of the figure represent the granulation flow 
field, and the thick curves represent magnetic field lines of flux tubes that are rooted in the 
intergranular lanes. The Ca II bright grains are thought to be located inside the flux tubes 
at heights of about 1 M m above the base of the photosphere. We suggest that the Ca II 
"straws" (lRuttenll2006l ) may be located at the boundaries between the flux tubes. 
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Fig. 3. — The temperature perturbation, AT, (about the initial state) at (a) 75 s, (b) 122 s, 
and (c) 153 s in a network element due to periodic horizontal motion at the lower boundary, 
with an amplitude of 750 m s~^, and a period of 24 s. The black curves denote the magnetic 
field lines, and the color scale shows the temperature perturbation. The white curves denote 
contours of constant /5 corresponding to /5 = 0.1 (upper curve), 1.0 (thick curve) and 10 
(lower curve). 
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Fig. 4. — The temperature perturbation, AT, (about the initial state) at (a) 75 s, (b) 122 s, 
and (c) 153 s in a network region consisting of 3 flux tubes. Wave excitation is due to 
periodic horizontal motion at the lower boundary, with an amplitude of 750 m s~^, and a 
period of 24 s. The black curves denote the magnetic field lines, and the color scale shows the 
temperature perturbation. The white curves denote contours of constant f3 corresponding 
to /3 = 0.1 (upper curve), 1.0 (thick curve) and 10 (lower curve). This figure is also available 
as an mpeg animation in the online version of the journal. 
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Fig. 5. — Field aligned (K) and transverse components of the velocity in a single tube, 
(panels a - b) and multiple tubes (panels c - d) at t = 120 s due to periodic horizontal motion 
at the lower boundary, with an amplitude of 750 m s~^, and a period of 24 s. The black curves 
are magnetic field lines, and the white curves are contours of constant /?, corresponding to 
(3 = 0.1, 1.0 and 10. In the weak-field region [B < 20 G) the decomposition of velocity into 
field-aligned and transverse components is ambiguous, and is not shown. 
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Fig. 6. — Close-up of the temperature perturbation AT and velocity field at time t = 101 s 
for (a) the central flux tube shown in Fig. HJ and (b) the single flux tube shown in Fig. [31 
The black curves are field lines, and the white curves are contours of plasma (3 [thick: (3=1] 
thin: /3 = 10). 
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Fig. 7. — Temperature perturbation, AT, (about the initial state) at (a) 76 s, (b) 101 s and 
(d) 140 s in a single tube due to a localised vertical periodic motion at z = 0, x = 2300 km 
(panels a- c) and x = 1900 km (panels d- f) with an amplitude of 750 m s~^, and a period 
of 24 s. The black curves denote the magnetic field lines, and the white curves are contours 
of constant /?, corresponding to /3 = 0.1, 0.3 and 1.0. 
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Fig. 8. — Field-aligned (V^) and transverse [Vn) components of the velocity at t = 140 s 
for a localized acoustic source at z = 0, x = 2.3 Mm (panels [a] and [b]) and x = 1.9 Mm 
(panels [c] and [d]). The black curves denote the magnetic field lines, and the white curves 
are contours of constant (3 {(3 = 0.1, 0.3 and 1.0). The velocity in the region corresponding 
to i? < 20 G is not shown. 



